Thermal conductivity of one-dimensional lattices with self-consistent heat baths: a 

heuristic derivation 
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We derive the thermal conductivities of one-dimensional harmonic and anharmonic lattices with 
self-consistent heat baths (BRV lattice) from the Single- Mode Relaxation Time (SMRT) approxima- 
tion. For harmonic lattice, we obtain the same result as previous works. However, our approach is 
heuristic and reveals phonon picture explicitly within the heat transport process. The results for har- 
monic and anharmonic lattices are compared with numerical calculations from Green-Kubo formula. 
The consistency between derivation and simulation strongly supports that effective (renormalized) 
phonons are energy carriers in anharmonic lattices although there exist some other excitations such 
as solitons and breathers. 
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Heat conduction exhibits diversified behaviors for one- 
dimensional lattices in terms of heat current j as the func- 
tion of lattice length JvQ, 0] • Harmonic lattice possesses 
ballistic heat transport [3[ in which the heat current is in- 
dependent of the lattice length, j ~ N°. Anharmonic lat- 
tices without external potential such as FPU-like lattices 
have anomalous heat conduction[4| in the sense j ~ N~ a , 
where < a < 1. Anharmonic lattices with external 
potential such as Frenkel-Kontorova (FK) and _6 A lat- 
tices show normal Fourier's heat conduction^, [g] since 
j ~ N . After intensive numerical simulations over the 
last decade, it has been commonly accepted that anhar- 
monicity and external potential are sufficient to induce 
Fourier's heat conduction law. In fact, Bolsterli, Rich 
and Visscher (BRV) have obtained in 1970 the Fourier's 
heat conduction by attaching external heat baths to ev- 
ery atom of a harmonic chain p}. In this so-called BRV 
lattice, the interaction between heat baths and harmonic 
chain resembles the external anharmonic forces. 

The BRV lattice plays a fundamental role in under- 
standing the Fourier's law in one-dimensional lattices as 
exemplified by recent studies [1, Q. In Ref.H, Bonetto, 
Lebowitz and Lukkarinen revisit the thermal transport of 
BRV lattice with a mathematically rigorous treatment. 
In addition, they also obtain the same result from the 
derivation of Green-Kubo formula. In Ref.0], they deal 
with the quantum BRV lattice and obtain the exact re- 
sult in the classical (high-temperature) limit. However, 
the phonon picture during thermal transport is not re- 
flected explicitly in their approaches. Harmonic lattice 
only allows phonons as the collective motion. There is 
no doubt that phonons should be the energy carriers for 
harmonic lattice. However for anharmonic lattices, there 
are some other collective motions such as solitons [lOj and 
breathers [llj existing as the consequence of anharmonic- 
ity. They are proposed as the energy carriers responsi- 
ble for some novel transport phenomena, e.g. thermal 



rectifying in asymmetric nanotubes[l2j, anomalous heat 
conduction in FPU-/3 lattice [l3| and so on , [l5| . An in- 
vestigation for anharmonic BRV lattices where external 
heat baths are attached to every atom of an anharmonic 
chain is desirable and this may shed some light on the 
role of energy carriers for anharmonic lattice. 

In this work, we study the thermal conductivities of 
BRV lattices from a different approach named as the 
Single-Mode Relaxation Time (SMRT) approximation 
from Boltzmann equation [lil Both harmonic and 

anharmonic BRV lattices have been considered. The re- 
sults from SMRT approximation are also compared with 
numerical simulations from Green-Kubo method. The 
harmonic BRV lattice is described as a one-dimensional 
harmonic chain of N atoms attached by independent 
Langevin heat baths. The Hamiltonian of this harmonic 
chain reads 
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where qi is the i'th atom's displacement from its equilib- 
rium position, u! 2 and 7 2 represent the strength of inter- 
atom and external potential respectively. For a chain 
with lattice constant a and atom mass m, we adopt the 
dimensionless units by setting a = m = 1. A periodic 
boundary condition q^+i = 9i is used. The heat bath 
attached to every atom is modeled by stochastic Langevin 
heat bath. The coupling strength is determined by the 
friction coefficient A in Langevin dynamics. Thus the 
equation of motion of the harmonic BRV lattice is ex- 
pressed as 

& = -o — (2) 

where £j is the Gaussian white noise with (£,i(t)) = and 
(0)> = 2Xk B TS(t). 
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According to the SMRT approximation [l6L 17 1, the 
thermal conductivity for one-dimensional system has a 
compact expression 



k = — I v (k)T(k)dk 
n Jo 
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where c is the specific heat, k the phonon wavevector, 
v(k) the group velocity, and r(fc) the phonon relaxation 
time of mode k. We focus on the classical transport pro- 
cesses where specific heat c = ks for harmonic lattice 
with Hamiltonian Eq. |T]). 

To obtain the group velocity v{k), we need to find 
the phonon dispersion relation. For Eq. |T]), one 
readily obtains the following dispersion relation uik = 

7 2 /w 2 . Thus 



w^/4sin 2 k/2 + v 2 for k G [0,tt] with v 2 
the phonon group velocity is obtained 
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If the system has no external potential, i.e. 7 = 0, the 
phonon group velocity is v(k) = to cos k/2. The existence 
of external potential shifts the phonon band off the zero 
point and reduces the phonon group velocity, as can be 
seen from Eq. ((4}. 

To analyze the phonon relaxation time r(/c), we con- 
sider the velocity auto-correlation function (qk(t)qk(®)) 
for mode k. Here qk is the counterpart of in mode space 
after some canonical transformation q^ = J"'., Sk%qi- S 
is the transformation matrix. We use symbol k to rep- 
resent wavevector and mode index simultaneously since 
it will not cause confusion. After canonical transforma- 
tion, the Hamiltonian can be decomposed into separated 
harmonic oscillators specified by a frequency Uk- Per- 
forming the transformation on both sides of Eq. ([2]), 
we obtain the separated equation of motion in mode 
space: q\ — ~io\q^ + £fc — \q^. After transformation, 
£fc still obeys the statistical properties of Gaussian white 
noise. Especially, the statistical properties of the fluc- 
tuations will not depend on the choice of the potential 
form. The above Langevin equation can be solved [ljjl 
by setting the first term on the right hand side to zero: 
q k (t) = q k (0)e- xt + f Q dt 1 e^-*') £, k {t'). Substituting the 
solution of <?fc(i) into the velocity auto-correlation func- 
tion, we have 



(q k (t)qk(0)) = (q 2 k (0))e 
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The ensemble average of the terms containing (t) van- 
ish due to the statistical property of Gaussian white 
noise. The velocity auto-correlation function decays ex- 
ponentially with a characteristic time of 1/A. Therefore 
we come to the conclusion that the phonons acquire a 
frequency-independent relaxation time r(k) = 1/A due 
to the coupling with the external Langevin heat bath. 




FIG. 1: (Color-online) Thermal conductivity k vs v. Solid 
curve is from Eq. (B), and bullets are numerical results 
from Green-Kubo method. The parameters are u — 1,A = 
0.25, T = 1 for a lattice with TV = 100. The Boltzmann con- 
stant ks is set to unity. 



As a result, the thermal conductivity can be obtained 
fromEq. © 
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Thus we have obtained the same exact result as in pre- 
vious works @, i, S from a heuristic derivation. The 
phonon picture for the thermal transport is explicitly re- 
vealed in our derivation. 

As we have mentioned in the introduction, it has been 
proved that the Green-Kubo formula renders the same 
result 0. Here we calculate the thermal conductivities of 
harmonic BRV lattice by Green-Kubo formula numeri- 
cally. The advantage of numerical Green-Kubo method 
lies in the fact that it can be easily extended to anhar- 
monic BRV lattice where no existing theoretical analysis 
has ever been applied. The thermal conductivity is cal- 
culated froml 81 
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(J(t)J(0))dt 



where J{t) — J2" =1 ji and ji = Qi ^-^^^ is the heat 
current within the chain. V(qi, qi+i) represents the inter- 
atom potential. The length and temperature indepen- 
dence of k has been verified by Green-Kubo simulations 
(not shown here). In Figfl] we plot the thermal con- 
ductivity as a function of v at lo = 1,A = 0.25, T = 1 



(7) 
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for a length N — 100. The numerical result for a 100- 
atoms chain has already approached the continuous-limit 
result. There is no surprise since the mean-free-path l(k) 
of phonons (l(k) — v(k)r(k) ~ oj/X ~ 4) is much shorter 
than the chain length N = 100. 

Next we shall consider the anharmonic BRV lattice. 
As we have seen above, the existence of external poten- 
tial only affects the group velocity. For simplicity, we 
only consider the anharmonic BRV lattice without ex- 
ternal potential term. First we consider the anharmonic 
case containing FPU-/3 chain. The Hamiltonian of FPU- 
[3 chain is 
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Due to the anharmonic interaction, the canonical trans- 
formation of Eq. ^ yields coupled equations of motion 
in mode space. In this case, we need to evoke the ef- 
fective phonon theory [3], which considers the phonon- 
phonon interaction as a kind of mean field. The anhar- 
monic Hamiltonian can be viewed as weakly coupled ef- 
fective phonons with a renormalized frequency ujk- Thus 
by ignoring the weakly coupled interactions within ef- 
fective phonons, we obtain the approximate equations 
of motion in mode space: q\ — —uj^qk + £fc — Ag^. 
The statistical properties of the fluctuations don't de- 
pend on the choice of the potential form. With sim- 
ilar analysis, the relaxation time for effective phonons 
is still the frequency-independent 1/A. The major an- 
harmonic effect comes from the group velocity of effec- 
tive phonons. According to the effective phonon the- 
ory, the effective phonon frequency ujk in FPU-/3 chain 
is proportional to the phonon frequency io k = 2sinfc/2 
with a temperature(anharmonicity)-dependent prefactor, 
i.e. Uk — s/acjk- The temperature-dependent co- 
efficient a has the following analytic expression a = 
l + / °° pe-PWSMWK 2 e-W 2 /2+0 4 /4) # with 
(3 = 1/fcgT. From this, we obtain the group velocity of 
effective phonons: v{k) = duik/dk = ^/a.cosk/2. 

The anharmonic interaction will also make the specific 
heat a function of temperature. This effect cannot be 
neglected when temperature is not very low. For a chain 
with N atoms, the specific heat is determined by c = 
(H) /NT. As a result of energy equipartition theory, the 
contribution to specific heat from kinetic energy is half 
of the Boltzmann constant The contribution from 

potential energy can be derived by using the equality 
NksT = J2iLi(li^fr)- Substituting the Hamiltonian, 
Eq. ([8]), into the above equality and using the definition 
of specific heat, one can obtain the specific heat for FPU- 

a 
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FIG. 2: (Color-online) Specific heat c of the FPU lattices as 
a function of temperature T. Solid curve is for FPU-/3 chain 
and dotted curve is for symmetric FPU-a chain. The high 
temperature limit values of specific heat, 3/4 for FPU-/? and 
5/6 for symmetric FPU-a chain (see Eq. (|lip h are marked 
by two arrows pointing to right. Boltzmann constant ks has 
been set to unity. 



The low and high temperature limit value of c are ks and 
§fe.B, respectively. The temperature-dependent c is plot- 
ted in Figf21 The monotonically decreasing specific heat 
as the function of temperature looks counter-intuitive as 
we all know that the specific heat of real material is a 
monotonically increasing function of temperature. We 
have to emphasize two things clearly here. First we only 
deal with classical specific heat. For harmonic chain, 
the classical specific heat is a temperature-independent 
constant ks- Second, we have introduced the anhar- 
monic interaction. The anharmonic interaction actually 
reduces the specific heat. Thus the "counter-intuitively" 
decreasing specific heat with temperature is a classically 
anharmonicity- induced result. 

The thermal conductivity of FPU-/? BRV lattice can 
be derived from Eq. (|3|) as 



c 

k = — 

7T 
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Thus, k depends on temperature via the product of c and 
a. To verify this prediction, we perform the Green-Kubo 
calculations of k. The results are plotted in FigO Solid 
curve is our prediction of Eq. (fTO]) and red solid circles 
are the numerical values of k. Our prediction agrees very 
well with the Green-Kubo calculations for a very wide 
temperature range of about four orders of magnitudes. 

It must be emphasized that our prediction is derived 
only by assuming effective phonons as the energy car- 
riers and ignoring the weak interactions between them. 
The consistency between theoretical analysis and numer- 
ical simulations clearly demonstrate that it is effective 
phonons, a kind of renormalized phonons due to anhar- 
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FIG. 3: (Color-online) Temperature-dependent thermal con- 
ductivity ft of anharmonic BRV lattices. Solid and dotted 
curves are from Eq. (|10[1 for FPU-/3 and symmetric FPU-a 
BRV lattices, respectively. Solid circles and hollow circles are 
the numerical Green-Kubo results for FPU-/3 and symmetric 
FPU-a BRV lattices, respectively. The horizontal line marks 
the value of k for harmonic BRV lattice. The simulations are 
performed at A = 0.25 for a lattice with N = 100. Boltzmann 
constant ks has been set to unity. 



monic interactions, carry heat energy in thermal trans- 
port processes for FPU-/? lattice. 

Besides FPU-/? BRV lattice, we also consider another 
simple anharmonic BRV lattice containing a symmetric 
FPU-a chain 
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The expression of specific heat 
obtained 



and coefficient a are 
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As we have discussed for FPU-/3 BRV lattice, the 
symmetric FPU-a BRV lattice also has a frequency- 
independent phonon relaxation time, i.e. r(k) = 1/A. 
The final expression of k is also the same as Eq. (|TU| . 
Here we should keep in mind that the specific heat c and 
a are now coming from Eq. (|12p . This prediction is plot- 
ted in Fig. [3] as the dotted curve. The numerical results 
(hollow circles in Fig. [3|) are in good agreement with this 
prediction. 

In summary, we have derived analytically the thermal 
conductivities of harmonic and anharmonic BRV lattices 



from the SMRT approximation. The derivation is heuris- 
tic and the phonon (effective phonon) picture is explicit 
during derivation. For harmonic BRV lattice, we obtain 
the same exact result as in previous works. For anhar- 
monic BRV lattices, we obtain the approximate results 
and compare them with numerical simulations from the 
Green-Kubo formula. The consistency between our the- 
oretical results and numerical simulations demonstrates 
that the effective (renormalized) phonons should be the 
fundamental energy carriers of anharmonic lattices. The 
contributions from solitons and breathers, if any, are neg- 
ligible at least for FPU-like lattices. 
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